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Abstract. We present an analysis of spectroscopic and photometric data for 
cool DQ white dwarfs based on improved model atmosphere calculations. In 
particular, we revise the a tmosp heric parameters of the trigonometric parallax 
sample of IBergeron et all d200lf) . and discuss the astrophysical implications on 
the temperature scale and mean mass, as well as the chemical evolution of these 
stars. We also analyze 40 new DQ stars discovered in the first data release of 
the Sloan Digital Sky Survey. 


1. Introduction 


Cool white dwarfs showing traces of carbon, either as neutral carbon lines or 
molecular C 2 Swan bands, are collectively known as DQ stars. Past analyses 
showed that these stars are helium rich with carbon abundances ranging from 

log C/He = —7 to —2 as determined from optical or ultraviolet. (IUE) spectro- _ 
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The presence of carbon has been successfully explained by a model in which 
carbon diffusing upward from th e core is brought to the photosphere by the 
deep helium convection zone ( Pelletier et al.lll986l h This model was shown to 
reproduce the observed carbon abundance distribution as a function of effective 
temperature fairly well for a helium layer thickness of logg(He) = Mn e /M* ~ 
—3.75, a value at odds with those obtained from evolutionary models of post- 
AGB stars that predict much thicker helium layers of the order of logg(He) ~ 
- 2 . 

More recently, iBergeron et al. (119971 1 and iBergeron et al. HI, hereafter 
collectively referred to as BLR) presented a photometric and spectroscopic anal¬ 
ysis of a large sample of white dwarfs aimed at improving our understanding of 
the chemical evolution of cool white dwarfs. Their sample includes several DQ 
stars that were analyzed with pure-helium model atmospheres, the only models 
available at that time. Recently. [Provencal et ali ( 2002 ) studied the importance 
of including carbon in the atmospheres of the DQ star Procyon B. Their effective 
temperature of T c g- = 7740 ± 50 K, derived from models including carbon, was 
sign ific ant l y coo ler than the earlier value of T e ff = 8688 ± 200 K obtained by 
[Provencal et al.1 (8 1997 1 on the basis of pure helium models. This large difference 
was attributed to an increase of the He - free-free opacity resulting from the 
additional free electrons provided by carbon. 
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The natural step following the comprehensive studies of BLR is thus to 
include explicitly metals and molecules in the model atmospheres. 


2. Theoretical Framework 

The m ode l atmosphere code is a modified version of that described at length in 
iBergeron et al.l ( 19951 1. which is appropriate for pure hydrogen and pure helium 
atmospheric compositions, as well as mixed hydrogen and helium compositions. 
In order to treat in detail white dwarfs of the DQ and DZ spectral types, we had 
to include metals and molecules in our equation-of-state and opacity calculations. 
Details will be provided elsewhere (Dufour et al. 2005, in preparation). Since in 
this analysis we are only concerned with DQ stars showing the C 2 Swan bands, 
we need only to include helium and carbon. Our model grid covers a range of 
T e g = 5000 to 12,000 K in steps of 500 K, log g = 7.5 to 9.5 in steps of 0.5 dex, 
and log C/He = —9.0 to —2.0 in steps of 0.5 dex. 


3. Analysis of the BLR Sample 

The first sample used for this study is drawn from the BLR analysis. It consists 
of 12 DQ stars with trigonometric parallax measurements, high signal-to-noise 
spectroscopy showing the carbon molecular bands, optical BVRI and infrared 
JHK photometry. The method used to fit the data is similar to that described at 
length in BLR, with the exception that in addition to the effective temperature 
and the solid angle n(R/D) 2 , we now have a third fitting parameter, the carbon 
abundance. Briefly, we transform the magnitudes at each bandpass into average 
fluxes. The resulting energy distributions are then fitted with the model fluxes, 
properly averaged over the filter bandpasses, using a nonlinear least-squares 
method. 

We begin by assuming a pure helium composition and log g = 8.0 to get 
initial estimates of T e g- and ir(R/D) 2 . The solid angle is combined with the 
distance D derived from the trigonometric parallax to obtain the stellar radius 
R. This radius is then converted into log g using an appropriate mass-radius 
relation for white dwarfs, and the procedure is repeated until a convergence in 
log g is reached. Next we turn to the spectroscopic observations and determine 
the carbon abundance by fitting the Swan bands at the values of T e g- and log g 
obtained from the fit to the energy distribution. This carbon abundance is then 
used to obtain new estimates of the atmospheric parameters from the energy 
distribution, and so forth. We repeat the procedure — typically five iterations — 
until T e g, log g, and the carbon abundance converge to a consistent photometric 
and spectroscopic solution. The logg valu es a re then converted into masses 
using the cooling models of lFontaine et al.l (2001)). 

Our spectroscopic fits for the BLR sample are displayed in Figure El the 
photometric fits differ very little from those already shown in BLR and will thus 
not be repeated here. 

We compare in Figure [3 the effective temperatures and masses derived from 
pure helium models (e.g. BLR) with those obtained here with models including 
carbon (see also Weidemann, these proceedings). We find that our revised ef¬ 
fective temperatures are reduced significantly with respect to the pure helium 
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Figure 1. Fits to the Swan bands plotted in order of decreasing tempera¬ 
tures from top to bottom. WD 0856+331 also shows atomic C I lines. 


solutions, in agreement with the conclusions of lProvencal et alJ 1 2002 ) for Pro- 
cyon B. The inclusion of carbon in the equation-of-state increases the number 
of free electrons and consequently the contribution of the He - free-free opacity. 
The resulting effects on the atmospheric structures can produce significant dif- 
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ferences in the derived atmospheric parameters. Since the T e g values have been 
greatly reduced, so are the model fluxes, and larger solid angles are required to 
fit the photometric data, which imply larger stellar radii or smaller masses, as 
shown in Figure [21 





Figure 2. Comparaison of effective temperatures and masses obtained from 
pure helium models (open circles ) and models including carbon (filled circles). 


An interesting consequence of our findings concerns the mass distribution 
of DQ stars. An intriguing question raised by BLR was the fact that the mean 
mass of DQ stars, (M) = 0.72 M 0 , was much higher than that obtained for their 
likely progenitors, the DB stars, bv iBeauchamn et alJ (1996), (M) = 0.59 M 0 , 
based on spectroscopic fits. Our improved analysis brings the mean mass of DQ 
stars down to ( M) = 0.61 Mg, in much better agreement with the spectroscopic 
mean mass of DB stars, and more in line with our understanding of white dwarf 
evolution. 


4. Analysis of the SDSS Sample 

Our second sample consists of 40 stars spectroscopically identified as DQ (hu¬ 
man ID) in the first data release of the Sloan Digital Sky Survey (SDSS). Details 
concerning these observations are described in iKleinman et all 1 2004 1 and ref¬ 
erences therein. Prior to the SDSS study, only 37 cool (T e g < 15, 000 K) white 
dwarfs with carbon features had been analyzed by various investigators, and 
the SDSS first data release alone has more than doubled the number of DQ 
stars known, allowing us to improve significantly the statistics of DQ stars. The 
method we use here to fit the SDSS data is similar to that described above, with 
the exception that the SDSS ugriz photometry is used instead. These photo¬ 
metric passbands cover the entire optical range from the UV atmospheric cutoff 
(3200 A) to the red sensitivity cutoff of the detector (~ 10,000 A). Furthermore, 











Analysis of DQ White Dwarfs 


5 


since trigonometric parallax measurements are not available for the SDSS stars, 
we assume a value of log g = 8.0 for all objects. 


5. Discussion 

Our results are summarized in Figure 0 where we show the carbon abundances 
as a function of T e ff for both the BLR and SDSS samples. We see that the 
carbon abundances form a relatively narrow sequence in the log C/He — T e g 
plane. The sequence defined by the SDSS sample overlaps the BLR sequence 
quite nicely, suggesting that the physical characteristics of these two samples 
are very similar. Since the BLR sample as a whole has a normal mean mass 
near 0.6 M 0 , the value of log g = 8 assumed for the SDSS sample is probably 
reasonable. It is also comforting to see that Procyon B (star sy mbol at 774 0 K), 
whh a precise mass determination of 0.602 M 0 determined by I Pi 1 oven cal et all 
( 2002 1. falls well within our sequence. Our results are thus entirely consistent 
with the idea that all DQ stars evolved from hotter DB white dwarfs. 



T eff (10 3 K) 


Figure 3. Observed carbon abundances as a function of T e s for both the 
BLR sample (filled circles ) and SDSS sample (open circles ); the star symbol 
shows the location of Procyon B. The error bars in the upper right corner 
indicate the average uncertainties of T e s (~ 150 K) and log C/He (~ 0.3 
dex). The dotted line represents our detection threshold for the C 2 Swan 
bands (or atomic C I lines for T e s > 10, 000 K). The solid and dashed curves 
correspond to the evolutionary sequences discussed in the text. 

Also shown in Figure El are the results from two different sets of evolution¬ 
ary model calculations, which predict the carbon abundances as a function of 
decreasing effective temperatures for various values of the thickness of the he- 
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lium layer g(He). The dashed curves represent the calculations of lPelletier et all 
( 1986f) at 0.6 Mq with, from top to bottom, logg(He) = —4.0, —3.5, and —3.0. 
As discussed in the Introduction, our abundance determinations would suggest 
values of logg(He) smaller than —3.5 for all DQ stars, much thinner than the 
values predicted by post-AGB evolutionary models. A new generation of evo¬ 
lutionary models computed by Fontaine & Brassard (these proceedings) yield 
quite different results, however, as shown in Figure 01 by the solid curves for 
models at 0.6 Mq with, from top to bottom, log (/(He) = —4.0, —3.0, and —2.0. 
With these new models, the observed carbon abundances in DQ stars are con¬ 
sistent with values of the helium layer thickness between logg(He) = —3 and 
—2, in much better agreement with the results from post-AGB models. 

Several objects in Figure 01 show carbon abundances that lie about 1 dex 
above the bulk of DQ stars. Among these, the only white dwarf with a measured 
trigonometric parallax is G47-18 (WD 0856+331; top filled circle) with an esti¬ 
mated mass of M = 1.05 Mq, by far the most massive DQ star in Figure 01 It 
is thus tantalizing to suggest that the DQ stars with larger-than-average carbon 
abunda nces rep resent the high mass component of the white dwarf mass distri¬ 
bution. iLiebert et alJ ( 2003) has suggested that the hot DQ stars discovered in 
the SDSS could all be massive, and could correspo nd to the missing high mass 
tail of the DB star mass distribution (fBeauchamp et al. ! 199(ii i. If both hypothe¬ 
ses are confirmed, the cool DQ stars with large carbon abundances observed here 
would simply represent the natural extension of the hot DQ stars identified by 
Liebert et al. Finally, the absence of neutral carbon lines in the optical spectra 
of DB stars with normal masses is naturally explained by our results. Indeed, 
the expected carbon abundances in hot (T e g > 10,000 K) helium-rich white 
dwarfs with logg , (He) between —3 and —2 fall below the detection threshold of 
the atomic C I lines according to Figure El 
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